NPS ARCHIVE 
1966 

ALVAREZ, F. 



FLUID CONTROL SYSTEM BEHAVIOR 
IN A RANDOMLY VIBRATING ENVIRONMENT 
by 

LCDR. FRANKLIN F. ALVAREZ, U.S.N. 

SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE 
DEGREE OF DOCTOR OF 
SCIENCE 
in 

MECHANICAL ENGINEERING 
at the 

MASSACHUSETTS INSTITUTE OF 
TECHNOLOGY 



AUGUST, 1966 



DUDLEY KNOX LIBRARY 
NAVAL ROSTORADUATE SCHOOL 
MONTEREY, CA t3f43-itt1 



1 



FLUID CONTROL SYSTEM BEHAVIOR 
IN A RANDOMLY VIBRATING ENVIRONMENT 
by 

LCDR. FRANKLIN F<, ALVAREZ, U.S.N. 

B.S., UNITED STATES NAVAL ACADEMY 
(1957) 

S,M., MASSACHUSETTS INSTITUTE OF TECHNOLOGY 

(1964) 

NAV.E., MASSACHUSETTS INSTITUTE OF TECHNOLOGY 

(1964) 

SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE 
DEGREE OF DOCTOR OF 
SCIENCE 



at the 

MASSACHUSETTS INSTITUTE OF 
TECHNOLOGY 



AUGUST, 1966 



7 (, 



DUDLEY KNOX LIBRARY 



- 2 - 
ABSTRACT 



FLUID CONTROL SYSTEM BEHAVIOR IN A RANDOMLY VIBRATING ENVIRONMENT 
LCDRo Franklin F, Alvarez, U.S.N. 

Submitted to the Department of Mechanical Engineering on August 15, 
1966 in partial fulfillment of the requirements for the degree of 
Doctor of Science, 



A fluid control system is analysed while operating in a 
vibrating environment. The introduction and transmission of noise 
signals in the way of unwanted pressures are examined for passive 
transmission line elements and an active flapper-nozzle control valve 
when subjected to environmental random or periodic translational 
vibrations. Methods are presented for predicting the levels and 
frequencies of induced pressures using distributed and lumped parameter 
techniques, including some damping effects. Analytic results evalu- 
ated with the help of a digital and analog computer are presented 
along with experimental data obtained for vibrating pneumatic trans- 
mission lines and for a small pneumatic control valve. Also included 
is a procedure for testing a given system’s susceptibility to vibra- 
tion induced noise sources. 

It is found that pressure peaks proportional to the fluid 
density, speed of sound in the fluid, and degree of excitation are 
generated in transmission lines over a wide band of frequencies. The 
peaks occur at the system resonant frequencies, and are due to induced 
fluid flow caused by pipe boundary nonuniformities (bends and area 
changes) moving relative to the fluid. The noise signals are excited 
by either longitudinal or transverse piping vibrations making the 
system susceptible to noise introduction regardless of its orienta- 
tion with respect to the vibrating field. 

Noise is introduced into the control valve by means of the 
vibrating supply line, inertial forcing of the valve flapper, and for 
high pressure systems, by vibrations of the transmission line to the 
load. The greatest noise transmission takes place at low frequencies 
and near the flapper resonance. Force feedback by means of a bellows 
arrangement decreases the amplitudes at the lower frequencies, but 
reinforces those near the flapper resonance and can cause noise power 
amplification. 
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1. INTRODUCTION 



hi backgroun d 

High performance hydraulic and pneumatic control systems are in 
wide use today because of the high speeds of response they offer. 

These systems can be found in missiles, aircraft, and ships, each of 
which is at times subject to random vibrations from either internal or 
external sources such as rocket combustion, engine vibration, or ran- 
dom seas. A fluid control system operating in this environment is 
subject to noise, unwanted pressure signals, excited by induced motion 
of the system structure and of the working fluid itself. How this 
noise is introduced into the system and how it effects the transmission 
of command signals are the two main questions this thesis will try to 
answer. Other important considerations are the dependence of noise 
generation on the orientation of the system to the vibration field, 
coupling between the fluid and structure, and transmission coupling 
into directions other than that of the primary vibrations. 

Studies of fluid control valve forces both steady-state and 
transient have been conducted in the past 2, ^ . However, none 

of the cases studied have included any environmental vibrations of the 
valve assembly itself with inertial forces being of significant impor- 
tance. 

Likewise, transmission line behavior has been investigated [a, 5, 6_ 
for a line that is in an inertial system. D’ Souza and Oldenburger 
included some coupling between the piping and fluid, but they still did 
not consider any vibrating environment induced motion. 
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Silberman [bJ did vibrate a tube with fluid flowing through it, 
but this investigation was solely to find the effect on pressure loss 
between two points, and it did not look at any pressures induced by 
the vibrations. 

In Lipner and Fay [ 9 J one finds the first mention of environmental 
induced pressures® This paper discusses the transient behavior of 
acoustic waves generated by a seismic shock. The simplified analysis 
begins with a pipe of infinite extent and ignores end conditions even 
when boundaries are established. Thus fluid resonances are not 
considered. 

Random vibration studies have been conducted investigating parts 
failures in electronics equipment pLO^ , structural response to random 
forcing [ll, 12, 13~] , and fatigue failure of materials [l2, 13, 14]] . 
Once again, however, there is a lack of information dealing with fluid 
power elements while operating in a randomly vibrating environment. 



1.2 OBJECTIVES AND PROCEDURE 

The objectives of this thesis are three fold. The first is to 
investigate the sources of noise production and its transmission in 
fluid systems subject to environmental vibrations. The second is to 
analyse the behavior of an active element, a control valve, when it is 
part of this system, and to be able to give procedures for estimating 
its susceptibility to this type of noise. The third is to be able to 
recommend the most profitable avenues for further investigation in this 



area. 
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In order to conduct meaningful experiments in conjunction with the 
analysis, the following procedures were adopted • Air was used for the 
working fluid because of the large quantities of pressurized air 
available and the simplicity of exhausting directly to the atmosphere. 

A nozzle-flapper control valve was decided upon as the active element 
due to its ease of construction and the availability of previous 
analytical information. Also, it was decided to use a choked orifice, 
a pure resistance, for the load as no further insight would be gained 
by having the transmission line terminate in a more complicated impedence. 
Both analog and digital computers were used to help interpret the 
analytic results. 

Throughout this thesis the units of pounds mass, pounds force, 
inches, and seconds are used for mass, force, length, and time respective- 
ly. This introduces in Newton's Second Law the factor g^, relating pounds 

2 

force to pounds mass, which is equal to 386 Ibm in per Ibf sec . 

The thesis contains in Chapter 2 a study of vibration induced 
noises in transmission lines subject to both longitudinal and transverse 
vibrations. The last section of the chapter contains a summary due to 
the relative importance of this chapter to the thesis. Chapter 3 
includes the analysis of a control valve operating in a vibrating 
environment, discussing the sources and transmission of noise. Also 
included is a procedure for testing a given system's susceptibility to 
vibration induced noises. The final conclusions and recommendations are 
presented in Chapters 4 and 5 respectively, and supplementary^ work is 



contained in the Appendix, 
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VIBRATION INDUCED TIUNSMISSION LINE NOIS E 

2.1 INTRODUCTION 

The basic system to be investigated consists of a supply line 
leading from a constant pressure reservoir to a flapper-nozzle 
control valve 0 The load controlled by the valve, in this case a pure 
resistance, is located at the end of a transmission line. Fig. 2.1. 

The entire system is assumed to be on a rigid base that can experience 
translational vibrations in all directions. Since the base is rigid, 
all supports, where joined to the base, experience the same motion 
independent of positions. The entire system is subject to induced 
vibrations due to the forced motion of the base that for this inves- 
tigation assumes the role of a rocket frame, plane fuselage, or ship 
hull excited by periodic or random accelerations. Due to the induced 
vibrations of the structure, noise is introduced into the system by 
means of motion of the control valve flapper relative to the nozzle 
and by the motion of the structure, containing the fluid, relative 
to the fluid. The valve flapper and actuator comprise a second order 
system and therefore are susceptible to noise when excited near the 
resonant frequency. Since their behavior is dependent upon the control 
valve properties, this source of induced noise will be investigated 
in Chapter 3. The rest of this Chapter will be devoted to fluid 
pressures induced by structure motion. 

2.2 CONCEPT 

When a vibrating pipe contains a fluid under pressure, the 
boundary condition at the interfaces is that the fluid has the same 
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motion as the piping. Thus for pipes experiencing vibrations along 
their longitudinal axis, if the contained fluid is of low viscosity 
such as air, the boundary layer is thin, and the main mass of fluid 
is at rest. However, the fluid at end junctions such as bends or 
area changes would be experiencing motion because of the longitudinal 
vibration of the solid boundaries. These junctions act as flow 
sources to the fluid in the same way as a piston oscillating in a 
stationary pipe. These flow sources excite fluid vibrations causing 
quite high pressures to be obtained at the resonant frequencies. 

For liquids, the boundary layer is thicker than for air, but the 
same general behavior takes place. Cavitation bubbles can form at 
the solid boundaries, however, reducing the effectiveness of the 
structure motion. 

When the induced vibrations are transverse, all the fluid is 
forced in a transverse motion along with the pipe. However, due to 
mechanical coupling the pipe sets up longitudinal vibrations at the 
bending resonant frequencies and at double the resonant frequencies. 
Then as before, any junctions act as flow sources to the fluid and 
impose longitudinal vibrations on top of the transverse motion. 

For this investigation, the base is assumed rigid and all support 
excitation is the same independent of position. However, in actual 
practice the piping structure can be excited independently at different 
locations especially if the distance between points is large. Thus 
each junction can act as an independent flow source, and the total 
fluid response will be the superposition of the responses to each of 
the sources. Likewise, if the excitation is both longitudinal and 
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transverse, the pressure response is the superposition of the two 
types, and each can therefore be treated separately in the analysis. 
These concepts will be expanded upon in the following sections. 



2.3 METHOD OF MODELLING VIBRATION INDUCED FLOW SOURCES 

The representation of the flow source concept by a physical 
system which can be modelled is shown in the following sketch. 








< > 



rs/ 




The representation allovjs a vibrating bend or area change to be 
modelled as two transmission lines forced at their junction by a 
flow source. Thus normal analytic techniques such as developed by 
Paynter [isj for energy transmission elements can be used. The 
forcing functions are then obtained when the additional conditions 
for the interaction between the fluid and piping are included in the 
analysis. These forcing terms which are peculiar to systems under the 
influence of environmental vibrations will be identified in a sample 
problem formulated in the next section. 
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2o4 DEVELOPMENT FO R LONGITUDINAL VIBRATIONS 

First to be considered is piping forced in a longitudinal 
motion. Since the fluid is forced at junctions it is the conditions 
at these junctions that are of most interest o Thus the analysis 
leads to the use of transmission matrices between junctions which 
can then be cascaded to obtain the overall effects. This method will 
be used here. A development of uniform, one-dimensional, linear, 
distributed system transmission matrices, along with extensions to 
systems in noninertial frames and with mean transport, is contained in 
Appendix A. 

In order to apply the transmission matrix, Eq , (A. 1-3), to the fluid 
in a transmission line it is first necessary to obtain a relation for 
the variables in differential equation form. That is. 



dx 



and 



dw 

dx 



= - Yp. 



This is necessary in order to identify the system parameters Z 
and Y when the power variables p, pressure, and w, mass flow rate, 
are used. 

Friction and damping will be ignored for the present as well as 
viscous transfer of motion from the pipe wall to the fluid. 

Using the continuity equation for a small elemental control volume 
of length dx in a uniform pipe gives, Fig. 2.2, 

“ ' ST 
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or 



0 



A 



ie.., 



f p 



dA 



( 2.4 - 1 ) 



ax ‘ at ‘ at 

where p is the fluid density and A is the flow area* 

When using the acoustic assumptions of small pressure fluctuations 
and an adiabatic process, the equation of state becomes 



dp = ■^ dp 

where 6 is the bulk modulus of the fluid. 

The change in area is found by conducting a stress analysis on 
a section of the thin-walled elastic pipe. This gives 




a(D - D,)dx = D.pdx 
0 1 1 ^ 

or 

da(D - D.) = D^dp 
o 1 i ^ 

! 

where a is the circumferential stress in the pipe and D and D, are 
the outside and inside diameters of the pipe respectively. 

Now the change in circumference with area is 

dC = |- dA, 
o 

and the stress-strain relation is 




where E is the modulus of elasticity and C is the circumference. 
Therefore the change of area with pressure becomes 
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dA = 



2A 



D 



dp . 



- » 

1 



Putting this and the equation of state into Eq • (2 •4-1) gives 



9w 

8x 



pA 9p 
"3 9t 



(2.4-2) 



where 



JL = i + 2 

3^ 3 E(D^/D.- 1) 

e 0 1 



(2 •4-3) 



It is to be noted that any compliance of the pipe is taken into 
account in Eq. (2.4-3) by decreasing the fluid equivalent bulk 
modulus and thus decreasing the speed of sound in the fluid. This 
effect is important in liquids which have a high bulk modulus. 
However, it can be neglected in air unless the pipe is very compliant 
and has a low modulus of elasticity. 

Now applying the momentum equation to the same elemental control 
volume gives, 



A ^ j 1 9vw 
- A — dx = — 



3x 



9 X 



dx + 



9w 

9t 



dx 



or 



9p _ V 9w ^ w 9v 1 9w 

9x g A 9x g A 9x g A 9t 
&0 



where v is the velocity of the fluid particles. 



From continuity. 



9w 

8x 



pA IJ + VA 



9p pA 9p 

9x B 9t * 

e 



(2.4-4) 
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Putting the equation of state in this expression gives, 



9v _ ^ 

9x 3 9t 3 3x * 
e e 



(2.4-5) 



Now putting Eqso (2.4-2) and (2.4-5) into Eq. (2.4-4) and noting 

from the equation of state that g — = g = c^, one can obtain 

^ o p ^o dp * 

2 



/ n iE. = + 

2 ^ 3x “ 2 3t g A 3t • 

c c ^o 



Since — << 1, this becomes 
c 



3p _ 1 3w 

3x g A 3t * 



(2.4-6) 



0 

Using the differential operator D = — , Eqs. (2,4-2) and (2,4-6) 

o t 

can now be written in the desired form. Thus, 



and 



- z« 

dx 



dx - ■ 



(2.4-7) 



where 



z » and Y = 

g A 2 

^o c 



Then s ^z/y', defined as the characteristic impedence, becomes 



~ g^A. 



and Y = propagation operator, becomes 

D 



Y = 



It is also desirable to express these relations in terms of a 
relative flow variable w^. Thus from Eqs. (2.4-7) one obtains 
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— ^ = - Zw = Zp ADx 
dx r o o 



and 



(2.4-8) 



dw^ 

dx 



^ - Yp, 



since w = w + p ADx . 

r o o 

x^ is the component of motion of the base in the direction of the 
longitudinal axis of the pipe, p^ is the mean fluid density, and 
Z, Y, Z^, and y bave the same values as above. 

The two sets of Equations, (2.4-7) and (2.4-8), can now be put 
directly into there corresponding transmission matrix forms, 



and 



3) and 


(A.; 


p(0) 




w(o) 




p(o) 




w^(o) 





coshyx 



sinhyx 

c 



Z^sinhyx 



coshyx 



p(x) 

w(x) 



Z ^sinhyx 



coshyx 

sinhyx coshyx 

•- c 



p(x) 



w^(x) 



(2.4-9) 



Z sinhyx (p ADx ) 
c o o 



(coshyx-1) (p^ADx^) 



(2.4-10) 

If transport is important, the transmission matrices can be 
obtained from Eqs. (2.4-9) and (2.4-10) by using the multiplicative 
factor as discussed in Appendix A. 

For the piping containing the fluid, the variables!, force, 
and absolute displacement will be used. Once more friction and 
damping will be ignored at present. Using Newton's Second Law for an 
element of pipe. Fig. 2.3, gives 

. M. 



Or 



9f _ 

~ - 2 



9t 



(2.4-11) 
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where A and p 
P 

respectively. 

becomes 



^ are the pipe cross sectional area and pipe density 
Now since f is positive when compressive, Hooke’s Law 



E 




(2.4-12) 



where E is the modulus of elasticity. 

Rewriting Eqs. (2.4-11) and (2.4-12) using differential operators, 
one obtains 



and 



dx 



- - 



ii 

dx 



- Yf 



where 



Z = 



p A 
P J 



and Y = 



1 

EA 

P 



(2.4-13) 



Thus, the character istic impedence and propagation operator are 

z 

CP P V g„ 



and 



Yp = D 






respectively. 



Once again it is desirable to obtain these relations in terms 
of a relative flow variable Thus Eqs. (2.4-13) become 

^ = - zc - Zx 
dx o 



and 



(2.4-14) 



dx 



= - Yf, 



» 





I 
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and Z, Y, Z , and v have the same values as above when the 
cp’ 'p 

substitution, 5 = S + x , is made. 

’ r o’ 

Then since Eqs. (2.4-13) and (2.4-14) are in the proper forms of 
Appendix A, they can be directly expressed as transmission matrices. 
Thus for the two cases, 

(2.4-15) 

Z sinhy x x 
cp 'p p 

(coshy x-l)x 
P P_ 

(2.4-16) 

Using the relations developed above it is now possible to 
investigate the behavior of fluid contained in a vibrating pipe. 
Consider first the simple problem of a pipe open to the atmosphere at 
one end and closed at the other. Let the open end experience vibra- 
tions in the direction of the pipe axis. 

r>yO 

L ■ ~i 

0 L 

The end conditions are, 

?j.(o) = 0, p(o) = 0, f(L) = -p(L)A, and w^(L) = p^AD5^(L). 



f (o) 




coshy X Z sinhy x 

P cp 'p 




f(x) 


5(0) 




siny X coshy x 

- cp ^ p J 




5(x) 



and 



f (o) 




CO shy X Z sinhy x 

P cp 'p 




f(x) 


^ 1 
o 

U 

^ I 




— sinhy x coshy x 

^ P P 

L cp ^ ^ J 




5j.(x) 



Using Eq. (2.4-16) gives 
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(1 - coshr )x 



5j.(L) = 



P o 



coshr 



f- sinhr p(L) 

cp ^ 

coshr 



where F L. 
P P 



Then from Eq. (2,4-10) one obtains, 

Z 



0 = fcoshrcoshr + (p A^D) sinhFsinhF H p(L) 
L p o pj 



+ Z sinhr (p AD) (1-coshr )x + Z sinhFcoshr (p ADx ) 
c o poc poo 



or 



p(L) = 



- Z p A sinh F Dx 
CO o 



coshFcoshr 4- — (p A^D)sinhFsinhr 
P ^ o 1 

cp 



where F = yL® Now if the excitation is measured in terms of a g of 



acceleration, which is 386 in per sec , this becomes 



p(L) = 



-386Z p A sinhFG 

C O 



D fcoshFcoshr + ~ p A^DsinhFsinhF "1 
L p Z o pJ 

cp 



The frequency response of 
for D, Thus, 



■P. (T-) 

G 



can be obtained by substituting io) 



The last two terms of this equation are the flow sources caused 
by the pipe’s vibrating. The next to last term is due to the 
closed end of the pipe moving relative to the open end, while the 
last term is due to the translational motion of the pipe. 
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pCL) ^ 
G 



- cp sin ~ 0 ) 
o c 



cos— OJCOS 
c 



.. P C A 

I p T o A 

Ig E g A 



, L . r p ^ 
— sin - oisinjl^^ Lu) 



o o p 

This equation shows that the pressure levels are proportional 
to fluid density, equivalent speed of sound, and excitation. High 
pressures will be induced at the frequencies at which the denominator 
rainimums occur. Relative minimums occur at the resonant frequencies 
of both the fluid and the pipe, and absolute rainimums occur at the 
overall system resonances. Damping would decrease the induced 
pressure levels somewhat, but would not change the overall appearance 
of the response. 

Another system might have the pipe mounted to the base by means 
of cantilevered supports. The pipe might then be supplied from a 
reservoir and exhaust through a choked orifice. Such a system. Fig. 

2M. 
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2.4, was analysed, and the theoretical frequency response for 
is plotted in Fig. 2.5. It can be seen that relatively high 
pressures have been induced by the vibrating environment at all the 
odd numbered fluid resonances. The large peak near 800 cps is due 
to resonance of the pipe on the cantilever supports. This also 
causes the fluid resonance near 850 cps to be much higher than the 
others as it is reinforced by the near resonant motion. 

When the cantilevers are stiffened by a factor of two, the 
resonant peak is moved far enough to the right that it is off the 
Figure. Also it is noticed that all but the first fluid resonant 
pressures have been decreased slightly. 



The Figure also contains a plot which includes a fluid mean 
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transport velocity of one tenth the speed of sound for the stiffened 
system. It is seen to have very little effect* This indicates that 
small values of transport are not too important. 

It is to be noted that in the Figure the resonant peaks are 
damped by the orifice resistance and by an assumed damping of 
approximately . 04 . 



2.5 DEVELOPMENT FOR TRANSVERSE VIBRATIONS 



The generation of pressures in pipes undergoing transverse 
vibrations is more complicated than in the previously developed 
longitudinal case. The pipe does exert a small component of force 
on the fluid in the direction of the pipe axis due to the bov/ed shape, 




- A |£-dx 

dX 




dx( 






but this does not produce any significant pressures as the slope 
(force component) is too small. Coriolis and radial accelerations 
are also insignificant for this particular configuration. Experi- 
mentally, however, relatively high induced pressures were observed 
along with longitudinal vibrations of the pipe ends even though they 
were braced. Significant end motion was observed at the pipe bending 
resonant frequencies and at double these frequencies indicating 
considerable coupling between transverse and longitudinal motion. 

One means of coupling can be attributed to the pipe maintaining 
approximately the same length while undergoing transverse vibrations. 
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The pipe curvature thus causes longitudinal motion of the ends in 

order to maintain this length. This type of coupling will be 

investigated now for a pipe under two different types of end conditions, 

pinned-pinned and c lamp ed-c lamped . In so doing the maximum transverse 

displacement of the pipe will be determined and used to find the 

change in projected length. The Bernoulli-Euler model for beams will 
* 

be used. 

Consider first a pipe pinned at both ends. Jacobsen and Ayre 
[it] give 



Y = c^cospx + c^sinyx + c^coshyx + c^sinhyx 

for the space part of the solution of the beam equation, 

.4 p A .2 ** 

Ljl + _R_e. iJL _ n 



(2.5-1) 



p A (1) 1/4 

where p = ( — ) , y(x,t) is the transverse displacement, 

SoEI 

Y(x) is the space dependent part of the displacement, and the c*s 
are the integration constants. I is the pipe cross-sectional moment 
of inertia. 



Reference [l^l P* 332 states th at i f the radius of gyration of the 
cross section of the beam, Y' "= VVA , is small compared to the 
wavelength investigated, then the Blrnoulli-Euler model is adequate. 

The normal beam equation without any fluid terms can be used as 
Long [is] has shown that the fluid has very little effect on 
transverse vibrations. This is especially true for air. 
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Now if the platform holding the end supports undergoes motion in 

a transverse direction, the end conditions will be 

d^Y 

at X = 0, Y = and ~ - 0, 

dx 

d^Y 

and at X = L, Y = and 0 

dx 

where y^ is the transverse motion of the supports. 

Putting these expressions into Eq. (2.5-1) and into its second deriva- 
tive gives 



and 




> 





tanh 



1!^ 

2 • 



Since this type of forcing causes resonances only at the 



symmetric, odd numbers, modes the maximum displacement will always 
occur at mid span» Thus writing Eq. (2.5-1) at x = y in terms of 
relative displacement, ^ becomes 



y.( T ) = 



2 cos 



yL 



2 cosh 



pL 



- 1 



2 “ 2 J 

Written in terms of g’s of acceleration the frequency response becomes 



>'r¥ 



386 



w 



2 cos — 



2 cosh 



pL 



- 1 



(2.5-2) 



This then gives the maximum relative transverse displacement of the 
pipe at various frequencies. 

Now following the same procedure for a pipe clamped at both 






m. 
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ends, reference [l7j gives 

Y = Cj^(cospx -f coshyx) ■+• c^Ccosyx - coshyx) 

+ Co(siniJx + sinhyx) -t- c, (sinpx - sinhyx) (2.5-3) 

3 4 

as a more convenient form for the space part of the beam equation 
solution when complicated end conditions are encountered. 

The end conditions are now 

at X = 0, y! = y and 4^-0, 
o dx * 

dY 

and at x - L, Y = y and - 0. 



Putting these into Eq. (2.5-3) and into its first derivative gives. 




C3 = 0, 



c 



2 



I 

y C 



o 2 



y 

o . sinuLsinhyiL + cospL - coshyL 
2 1 - cospLcoshyL 






and 



^4 = 



y c, = 
•'o 4 



sinyiLd - coshyL) 4- sinhpL(l - cospL) 
1 - cosyLcoshpL 



As before, writing Eq . (2.5-3) at x = ~ in terms of relative 
displacement and g’s, gives for the frequency response 



Yr# 



386 

0)2 



l/2(cos + cosh 






Caicos ^ - 



cosh 



yL. 

2' 



+ (sin ^ - sinh ^) - ij 

Reference [l7^ states that the shape of a pinned-pinned beam forced 
at the supports at resonance approaches the normal-mode shape, 
sin - , corresponding to the resonant frequency. Since it is 



resonance that is of interest the pinned-pinned pipe displacement is 
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taken to be 



, HTTX , 

y = y sin — — sinu)t 
•'r •'ro L 



(2,5-5) 



where n = Is the number of the mode and y is the relative 

9 9 9 

maximum displacement obtained from Eq. (2, 5-2) for a given assumed 
damping . 

Using this equation the change in length of pipe with pinned- 
pinned end conditions can be found as follows: 

4^ 




or 



dS = 



1 2 
I ^ ± / — L 

2 ^ 8x ^ 



dx 



as 



9y 



8x 



- << 1 , 



Then 



/ 

L r- 



S = 



1 2 
1*2 < 



dx . 



Now the assumption of this section is S = L so that 



/ 



2 

1 * 1/1 



dx 



However since Lf^L , with little loss in accuracy this can be written 
L ^ 



L = 



2 

1 * 1/1 > 



dx 



or 






Is 
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2 



(2.5-6) 



Putting in Eq. (2.5-5) gives 
3y 

•^r mr mrx 

^ * ^ro T T 



and 



9 V 

, ^ X f mr,.2 2 niix . 2 ^ 

^ IT ^ T ^ IT 



— (y niT 2 . 2 . . 2n7TX ^ 

2 -^ro — ) sin 03t(l + cos — ) 



Therefore, 

AL 



1 /^ro^^ ^2 ,2 f . 2n7Tx ^ , 

( — - — ) sin 6)t V (1 + cos — — ) dx 



1 ,2 . 2 . ^ L . 2mrx . 

4 ^~~L~ ^ ^ T~ ^ 



or 



AL = L ( ) sin cot* 



(2.5-7) 



Now following the same procedure for the d amped-clamped pipe 
the normal-mode shape jj.7^ is 

Yr - jjsinhpx - sinpx + a(coshyx - cospx) 



sincot 



where 



a = 



= sinhpL - sinpL 
cosyL - cosyL 



and the maximum relative displacement, is obtained from Eq, 

(2.5-4). 



Differentiating, squaring, and integrating the expression in Eq . 
(2.5-6) gives 
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pL 



ro .2, 



AL = ^ ) L ^l/2(l/2sinh2pL + pL)- sinpL cospL + coshpLsinpL) 

+ l/2(yL + sinpLcosyL)+ 2a |^1/^ cosh2pL + 1/2 + 1/4 cos2yL 
+ l/2(sinhyLsinyL - coshyLcosyL) - 1/2 (coshyLcosyL+sinhyLslnyL) 






l/2(l/2 sinh2yL -- yL)4*(coshyLsinyL - sinhyLcosyL) 



+ 1/2 (yL - sinyLcosyL) 



>]} 



sin u)t. 



( 2 . 5 - 8 ) 



Since Eq. (2.5-8) only holds for the resonant frequencies which 



are 

yL)^ = 4.730, 
yL)^ = 10.996, 
yL)^ = 17.279, 
and yL)^ = 23.562, 

and since cosh 10.0 = sinh 10.0 = 11,013.2 so that a = -1, the Equation 
can be simplified for all the resonant modes above the first. Cancell- 
ing terms it is found that 

- 1 - 1/2 cos2yLj sin^wt. (2.5-9) 

Thus the change in pipe length is determined for the two types of 

end conditions, and Table 2.1 can be constructed using Eqs . (2.5-7), 

(2.5-8), and (2.5-9). 

2 

Now since sin cjt = 1/2 - 1/2 cos2u)t, the experimentally found 

end double resonant frequencies can be attributed to this means of 

coupling, and the end accelerations can now be calculated as soon as 

the pipe length and ratio of maximum displacement to length are known. 

^ro 

Some typical values of were obtained using Eqs. (2.5-2) 

and (2.5-4) for a sinusoidal platform acceleration of one g and an 
approximate damping ratio of .04, and are presented in Table 2.2. 



AL = 



yL 



ro 



l^pL 
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TABLE 2.1 

Pipe Longitudinal End Movement Due to Bending 



Mode 

1 

3 

5 

7 



P inned -P inned 



AL 



2.47 ( 




)^L 



. 2 

sin (jJt 



22.21 

61.68 

120.90 



Mode 

1 

3 

5 

7 



C lamp ed -C lamp ed 



, . ^ro .2. .2 

6,35 ( — ) L sin mt 

57.71 



144.95 

271.69 



TABLE 2.2 

Ratio of Mid-Pipe Transverse Displacement to Pipe Length 
(pipe is 1/4" copper tubing filled with air at 200 psia) 



P inned -P inned Clamped-Clamped 



Mode 


Cps 


L 


!ro 

L 




Mode 


Cps 


L 


^ro 

L 


1 


34 


23 


7.5 X 


10-3 


1 


76 


23 


1.9 X 10"3 


3 


218 


27 


4.5 X 


10"^ 


3 


296 


27 


1.1 X 10 ^ 


5 


606 


27 


2 .3 X 


10"3 


5 


733 


27 


0.6 X 10"3 


7 


240 


60 


6.2 X 


10"^ 


7 


276 


60 


1.6 X 10 ^ 




Now u 


sing the 


results 


of Tables 2,1 


and 2 . 2 


some 


typical 



values of the ratio of end double frequency accelerations to base 
accelerations are calculated and given in Table 2.3. 
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TABLE 2.3 

Ratio of Pipe End Longitudinal 

Acceleration to Platform Transverse Acceleration 
(pipe is 1/4** copper tubing filled with air at 200 psia) 





Pinned- 


Pinned 




Mode 


2xCps 


Accel 


^end 

^base 


1 


68 


.76 




3 


436 


1.18 




5 


1212 


.07 




7 


480 


.0033 





Clamped~C lamped 



Mode 


2xCps 


Accel 


1 


152 


®base 

.62 


3 


592 


.34 


5 


1466 


.016 


7 


552 


.00065 



The values from the Table will be a little high due to the 
restraining action of the end supports® However, they show that the 
end accelerations are of a relatively high value for the first and 
third modes, but that they drop off quickly after that. This trend 
and the above values will be compared with experimental results in 
Section 2.8. 

While the above analysis gives valid results for the twice 
resonant frequency end motion, only a qualitative explanation can be 
given for the even stronger coupling at the resonant frequency. 

This effect can be due to nonsymmetry in the system since 
an unsymmetrically supported pipe can be forced into longitudinal 
motion by transverse vibrations as shown. 
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Thus any nonsymmetry could cause the longitudinal vibrations 
to begin, and once started, they very probably would be maintained. 

If the platform is excited with components along both the 
longitudinal and transverse axis, the longitudinal end motion will 
be the additive effect due to botho 

2.6 DAMPING AND ITS EFFECTS 

The previous sections have shown that all the important pheno- 
menon occur at resonances. Therefore, in order to have any feeling 
for the amplitudes involved it is essential to know the amount of 
damping. So far very little has been said about this for several 
reasons. First, most of the mathematical developments become con- 
siderably more complicated when losses are taken into account. 
Secondly, it is difficult to determine what form the losses should 
take, i.e. internal friction, external friction, etc. Lastly, the 
correct values for the parameters are very Ward to determine, and 
sometimes vary by as much as a factor of twenty. An example would 
be support damping of beams which changes greatly depending on the 
tightness of the connection. Thus the most satisfactory way to gain 
an understanding of the damping involved in a given situation is 
experimentally . 
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* 

The definition for lightly damped systems, 
to 

A(o = 2 Ccl) = , (2.6-1) 

n Q 

where A(o is the bandwidth, C is the damping ratio, and u)^ is the 
resonant frequency, gives a means of experimentally determining damp- 
ing by measuring the bandwidths of resonant peaks . This measurement 
is made at the half power point which is .707 of the peak amplitude. 

A plot of Eqo (2.6-1), Figo 2.6, shows that as damping increases 
the peak amplitude decreases and the peak broadens out . Brown [l93 
had also shown theoretically that fluid damping is frequency dependent 
and increases significantly as (o goes up. 

Eq. (2.6-1) gives good approximate values of damping from 
experimentally obtained amplitude verses frequency plots. However, 
it is desirable to be able to predict the effects of damping on 
amplitude in analytical discussions as well. This can also be done 
by using Eq. (2.6-1) and experimentally obtained values of damping. 

The zero damped infinite peaks of a plot can be rounded off by making 
its shape conform to the correct shape, in Fig. 2.6, for the given 
damping. This will give a good estimation of the peak amplitudes 
if the values of damping used are close to that actually experienced 
by the system. Good repeatable damping values for both the fluid 
and mechanical systems were obtained experimentally as discussed in 
Section 2.8. 



The fluid and mechanical systems in this investigation are all 
lightly damped . 



*• 



'* 9 

4 luie 







I 
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2,7 RANDOM EXCITATION 



The developments in sections 2.4 and 2.5 were conducted using 
sinusoidal excitation to the base. Although rockets, airplanes, 
and ships are often subject to narrow-band periodic disturbances, 
a majority of the time the excitation will be of a wide-band random 
nature. A convenient means of characterizing a wide-band random 
signal is to use a stationary process which has a uniform, band- 
limited spectral density. This is known as white noise. 




White noise is a close approximation to the desired physical process, 
and as long as the frequency band of interest of the system to be 
investigated falls between and it will give the desired results. 

One can develop the relation for a system response spectral 
density when the excitation spectral density is known. A development 
is contained in Crandall and Mark [l4^ , pages 70-71, and the 
transformation notation used is 



S((d) 



2tt 




R(i)e 



-iu3i 



dx 



(2.7-1) 



where S(w) is the spectral density and R(t) is the autocorrelation 
function. The results given are 



S (u3) = H(-a))H(ai)S (u)) 

y X 
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or 

2 

S (u) = |h(w)| S (o)), (2,7-2) 

y ^ 

where S^(co) is the response spectral density, S^(o)) is the exc'ftation 
spectral density, and H(o)) is the frequency response of the system. 

The response spectral density is relatively simple to obtain from 
Eq. (2 •7-2) since the frequency response for the system can be found 
using the techniques developed in Sections 2o4 and 2.5, and S^(oa) is 
si.mply 

s 0). < io)! < 0)^ 

o 1 — ' ‘ — 2 

0 elsewhere 

for the excitation assumed in this investigation. 

It is to be noted that the spectral density covers both positive 
and negative values of frequency and is in the units of quantity 
squared divided by circular frequency o Therefore if experimental 
work is being done and only positive values of frequency are used 
and one is working with cycles instead of radians, then in order for 
the excitations to have the same value a factor of 4n must be used, or 

W (f) = 4TrS^(o)). 

Then 

Wy(f) = |H(f) (2.7-3) 

where and are the excitation and response spectral densities in 
terms of cps . 

2.8 EXPERIMENTAL RESULTS 

In order to verify the theory developed in the previous sections. 







m 
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experiments were conducted by actually vibrating pipes. The pipes 
were vibrated longitudinally and transversely with both periodic and 
random excitation while containing air under pressures of 100 and 
200 psia. In some cases the air had a small mean transport and in 
others it had no flow. Some preliminary work was done with water in 
the pipes, but this was discontinued due to the difficulties 
encountered in trying to remove the entrapped air bubbles. The 
piping used was standard one quarter inch copper tubing. The vibra- 
tion generator was built from a high fidelity speaker, and the 
pressures and accelerations were obtained by using piezoelectric 
transducers. Data was taken up to 1000 cps since above this fre- 
quency the system response is almost nil. 

Figs. 2.7 and 2.8 show the piping arrangement for longitudinal 
and transverse vibrations respectively. One end of the piping is 
connected by flexible hose to a constant pressure reservoir while the 
other end terminates at the pressure transducer. Fig. 2.9. Additional 
details of the experimental procedures and instrumentation are 
contained in Appendix C. 

Fig. 2.10 shows good agreement between the experimental and 
theoretical frequency responses for a system of total length 51.6 
inches forced by longitudinal vibrations while containing air under 
200 psia pressure. There is a choked orifice near the transducer 
end of .032 inches diameter which gives a mean velocity of 14 feet 
per second or about one one-hundredth the speed of sound. The Figure 
shows fluid pressure peaks at all the odd resonances due to a cosine 
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function characteristic equation, and a peak near 850 cps due to the 
longitudinal resonance of the pipe<. The resonances are not multiples 
of each other because of the frequency dependence of boundary layer 
friction which causes the boundary layer to be thinner at higher 
frequencies. The analysis used in this investigation did not take 
this change in boundary layer thickness into account so that a 
frequency correction to the equivalent speed of sound must be made 
from Figo 2. 11, It is necessary to use this Figure to obtain the 
correct value of c in order to compute the correct values of the 
resonant frequencies , 

From Figo 2^10 it is seen that high enou^ pressures are 
generated per g of acceleration at the first several resonances to 
be of concern if the system were subjected to any substantial 
environmental excitation. Also, from the first two peaks, using 
Eqo (2o6-l), the fluid damping ratios are found to be approximately 
.03 and .04 respectively. 

Fig. 2ol2 is the frequency response for the pipe longitudinal 
accelerations at the nonforced end. This Figure shows good agreement 
in resonant frequency between the theory and experimental results, 
but the agreement in bandwidth is not too good due to limitations in 
the equipment resolution capabilities » The resolution at high 
frequencies is only one tenth as good as at low frequencies. The pipe 
end had to be mass loaded in order to bring the first resonance down 

* 

This Figure was taken from Brown 19 , where c is the actual speed 
of sound, c is the isentropic speed of sound, a is the pipe radius 
in inches, and is the kinematic viscosity in in^/sec. 




4 
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to below 1000 cps. It is seen that relatively high accelerations can 
be obtained which if they occur ed at low enough frequencies or if 
they corresponded to fluid resonant frequencies would cause pressure 
signals large enough to be of great concern , From this Figure the 
pipe damping ratio can be seen to be about »04o 

Now the same system was excited with white noise input. Figs. 

2.13 and 2.14 show the pressure response for reservoir pressures of 
100 and 200 psia respectively. Notice that the two Figures are 
exactly the same except the first is everywhere between five and 
seven db less than the second. This is exactly what would be expected 
as the induced pressure is proportional to density which for a gas is 
proportional to pressure. Thus when the pressure is doubled there 
should be a rise of six db in induced pressure noise, Fig. 2.15. 

These two Figures also show the fluid resonant frequencies at 60, 

200, 355, 455, 600, and 740 cps. Notice that the resonances at 455 
and 740 cps are much smaller than the others. This showed up in the 
theory also. The resonance at 600 cps is higher than it would 
normally be since although white noise was put into the vibration 
generator, the generator had resonances around 370 and 530 cps which 
show up throughout the graphs. The peak at 850 cps is due to the pipe 
longitudinal resonance. 

The several runs attempted with water filled pipes produced 

pressure levels approximately four times as high as for the 200 psia 
- 



An explanation on reading db graphs is given in Fig. 2.15. 
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air for a given excitation. The increased density and speed of 
sound would predict level increases by a factor of ten. However, 
increased damping lowers the resonant peaks so that the factor of 
four seems reasonable with the theory* 

There is a misleading upward trend to all the graphs made by the 
analyser and recorder* This is because the filter bandwidth is 
always a constant percentage of the mid frequency being swept, in all 
these cases 6%. Thus the filter bandwidth doubles each octave, and in 
the case of white noise, twice as much signal is being averaged* 

g^ 

”oc t 

Thus there is a 3 db rise per octave since 10 log — — = 10 log 2 = 3. 

g 

This 3 db rise can be seen in Fig. 2.16 which is the random white 
noise input to the vibration generator. 

Fig. 2ol7 shows the acceleration at the nonforced end of the 
pipe* The peak near 900 cps is due to the pipe longitudinal resonance, 
and the two peaks at 370 and 530 cps are the resonances of the shaker. 

In order to better understand these last several results. Fig. 2.14 
has been replotted in Fig. 2.18. In this Figure the 3 db rise has 
been eliminated and the pressure has been normalized by the input 
acceleration. Theoretical results are also plotted and agree 
favorably* This Figure shows the big drop in pressure after the 
first two resonances, and also shows that the resonances at 455 and 
740 cps are almost nil. The resonance at 600 cps is much higher than 
predicted as explained before. Because of the wideness of the filter 
at higher frequencies, the resonance at 875 cps is masked in the pipe 
resonance in the experimental results. 
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Before leaving longitudinal excitation, Fig. 2.19 is included to 
show that the fluid resonances are dependent on the total pipe length 
and not on the forcing flow source positions. For this case, the 
distance between flow sources (pipe bends) was left the same as above 
while increasing the total distance between the reservoir and pressure 
transducer. It can be seen that the fluid resonances are at lower 
frequencies, 40, 156, 250, 350, 450, and 550 cps, corresponding to the 
greater length. As before the high peak at 550 cps is due to the 
shaker resonance and the peak near 850 cps is once again the pipe 
resonance . 

Figs. 2.20 and 2.21 show the pressure, response to random transverse 
excitation for c lamp ed-c lamped pipes of lengths 27 and 23 inches 
respectively. The overall lengths from reservoir to pressure transducer 
are 46.1 and 42.1 inches, and in both cases there is zero mean flow. 
These pipes were not excited by inertial forces as the vibration 
generator was not large enough to move the whole system. Instead the 
pipes were forced directly at a point near one end. This method 
proved satisfactory. However, since the forcing was not S3mimetric as 
in the case of inertial loading, all the resonances are present and 
not just the odd numbered ones. 

Notice the difference between the two Figures. The first has many 
small resonances while the second has several large ones. This is 
because in the case of the 23 inch pipe, the fluid resonances and pipe 
bending resonances coincided. The two graphs are similar in that the 
first resonances are large while the higher ones fall off quickly. 

This is more noticable if the 3 db rise per octave is subtracted. This 
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fall off was predicted and is much more severe than for longitudinal 
excitationo 

Now looking at Fig. 2.20 in detail, the first peak at 44 cps is 

S'? 

excited by the first bending resonance for a pipe of this length with 
clamped-pinned end conditions. Therefore, even though the pipe ends 
where clamped, forcing the pipe near one end excited this mode of 
bending also since the support moved torsionally. The peak at about 
86 cps is due to the double resonant frequency as predicted. The 
second and third resonant frequencies cause the peaks at 130 and 
260 cps respectively, but both their double frequencies are masked 
by other peaks. The one is masked by the third resonant frequency 
and the other by the shaker resonance around 530 cps. The peak at 
57 cps IS excited by the first resonant frequency for claraped-clamped 
end conditions. Its double frequency shows up near 120 cps. The 
fluid resonances are at 68, 220, 368, 515, 660, and 803 cps respectively. 
Only one, 517 cps, corresponds with any of the acceleration peaks. 

In Fig. 2.21 the fluid resonances can be observed at 80, 242, 

410, 585, 741, and 920 cps. The bending resonances for clamped- 
clamped end conditions are at 73, 216, 420, and 691 cps. Since the 
first several correspond almost exactly with the fltild resonances, the 
peaks are seen to be relatively large. A small peak is seen at the 
first double resonant frequency also, about 152 cps. There is some 
motion of this beam in the clamped-pinned mode, but since these 
* ^ 

The pipe bending resonances were computed by directly using the 
formulas listed in Den Hartog . 
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frequencies do not correspond to any fluid resonant frequencies, the 
small peaks are masked out. 

Fig* 2.22 shows the longitudinal accelerations at one support of 
this pipe* For the clamped-clamped mode the first resonance and its 
double frequency show at 73 and 158 cps respectively. Likewise for 
the claraped-pinned mode the first two resonances and thpir double 
frequencies are seen at 58, 111, 184, and 360 cps* The peak at 39 cps 
must be due to the first resonance of the pinned-pinned mode. 

In these last three Figures it is seen that the primary bending 
frequencies do excite pressures through longitudinal motions. Also, 
the double bending frequency accelerations are present, but they are 
of smaller magnitude and only show up for the first few resonances. 
Table 2*4 compares some of the experimental and theoretical results. 

TABLE 2.4 

Comparison of Experimental and Theoretical Ratios of Pipe End 
Longitudinal Acceleration to Platform Transverse Acceleration* 

(pipe is 1/4^* copper tubing filled with air at 200 psia) 





Clamped-C lamped 






Clamped-Pinned 




Cps 


F ®end 

^ 

®base 


Theo 

®base 


Q2S. 


®end 

^ e 

®base 


Theo 


120 


.32 


.74 


98 


.40 


.74 


156 


.28 


.62 


111 


.23 


.71 



^end 

7 

^base 



The experimental data had to be converted to g’s per base acceleration 
since the pipe was not inertially forced. This was done by using 
Eq. (2.5-2) and a similar one for clamped-pinned end conditions. 
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This Table shows the theoretical results to be a little high 
because the end support stiffness in reality places some restraint 
on the motion. However, the results are of the right order of 
magnitude and show that the techniques of Section 2.5 could be used 
to determine an upper bound for the end accelerations. 

On the other hand the magnitude of the longitudinal end accelera- 
tions due to the primary bending resonances can not be determined at 
the present time. Considerably more experiments with various end 
conditions and base supports will have to be conducted in order to 
be able to predict this effect. The most that can be said now is that 
the first resonance is the most important and seems to induce a 
pressure peak of approximately the same magnitude as that induced by 
comparable longitudinal excitation. 

The last two Figures included are to show the combined effects 
of both longitudinal and transverse accelerations. The pipe was 
excited near one end at an angle of 45 degrees to its axis, giving 
both types of forcing. Fig. 2.23 shows the pressure response to have 
the general upward trend, on the recorder graph, of longitudinal 
excitation, with fluid resonances near 80, 242, and 400 cps. However, 
the bending resonances are also present. For the clamped -pinned mode 
the first resonance and its double frequency are at 59 and 121 cps 
respectively, but those for the clamped-claraped mode are masked behind 
peaks. They can be seen, however, in Fig. 2.24 at 77 and 154 cps 
along with the second bending resonance at 216 cps. These two Figures 
show that both longitudinal and transverse vibrations can excite fluid 
pressures at the same time, making the system susceptible to noise 



ilTaffahiri < 
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regardless of the orientation of the system with respect to the 
vibration field. 

Finally, Table 2.5 is included with a summary of the experimentally 
determined damping ratios. 



TABLE 2.5 

Experimental Damping Ratios 



Fluid-no 


1 mean 


flow Fluid-mean flow 


Pipe-long 


Pipe- 


trans 


(200 


psia) 


(200 psia. 


14 ft/sec) 


d/4" copper 


tubing, air 


at 200 psia) 


Cps 




Q2S. 


£ 


C£S £ 


Cps 


£ 


45 


.03 


42 


.04 


850 ,04 


73 


.04 


80 


.03 


155 


.04 


900 ,04* 


78 


.05 


232 


.04 


197 


,04 








390 


,04 


350 


.04 









It can be seen that the values are fairly consistant, having 
a slight increase with frequency. Also, the fluid damping with no 
flow is smaller than with flow since the orifice adds resistance when 
it is in the system. 



2.9 SUMMARY 



In Chapter 2 it has been shown that piping systems containing 
fluid can excite large pressures v/hen their supporting environment is 
subjected to vibrations. Energy transfer to the fluid is by means of 



i: 

This value was corrected from a recorder output value of .07. 

Since the filter bandwidth is approximately 50 cps at this point 
it gives a misleading broadening of the peak which must be corrected. 








i 



- 46 



flow sources caused by pipe boundary nonunif ormities such as bends 
and area changes moving relative to the fluid. When subject to a 
white noise environment the pressure spectrum is composed of numerous 
resonant peaks over a wide band of frequencies, the peaks being due 
to the various resonances in the system. 

The piping supports can be forced longitudinally, transversely, 
or in combination. In all cases it is the longitudinal motion of the 
end nonuniformities that forces the fluid since the pipe transverse 
bending vibrations also cause longitudinal motion of its ends. These 
longitudinal end motions are excited at the bending resonant frequencies 
and at twice the bending resonant frequencies. 

The induced pressures due to longitudinal environmental vibrations 
were predicted and the experimental data gave excellent agreement with 
the theory. The excitation coupling to the twice bending resonant 
frequency end motion was first observed and then a theory developed 
which allows one to predict pressure levels with a fair degree of 
accuracy. The coupling effect at the primary bending resonant 
frequencies was not anticipated, but the theory for longitudinal 
vibrations can be applied once the coupled longitudinal motion has 
been measured. The experimental results show this coupling to be such 
as to induce pressures of the same order of magnitude as those induced 
by comparable longitudinal excitation. 

Since both longitudinal and transverse accelerations of the 
piping supports induce fluid pressures, the system is susceptible to 
noise regardless of its orientation to the vibration field. Motion 
components in more than one of the directions simply cause an additive 







1 



- 47 



effect to the longitudinal forcing motion. 

The level of pressure induced for a given excitation was seen 
to be proportional to the effective speed of sound in the fluid and 
to the density which is proportional to the pressure for compressible 
fluids. For the particular system investigated, air at 200 psia, 
pressures of approximately .6 psi per g of environmental acceleration 
were observed for the first several resonant peaks. These pressures 
would have been even greater for fluids of higher density as was 
indicated by the scanty data obtained for water filled pipes. 

Reference [iB^ contains acceleration spectra for aircraft, rockets, 
and ships which can be used to estimate the levels of induced pressures 
that systems in actual operation can experience. 

Extremely high pressures can be excited when a fluid resonance 
coincides with that of the piping system. For a more complicated 
system it is impossible to distinguish between these two types of 
resonances, but there are still absolute system resonant peaks that 
are higher than Che others. This phenomenon can cause significant 
induced pressures at the higher frequencies even though the fluid 
resonant peak levels normally fall off sharply after the first two 
or three. 

Experimental results show the fluid resonant peaks to be fairly 
sharp, and both the fluid and piping damping ratios appear to be .04 
for the cases tested. 

The good agreement between the theoretical and experimental results 
indicates that the assumption made neglecting the viscous boundary 
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layer transfer from pipe to air is permissible over a wide band of 
frequencies . 

Reducing the production and transmission of these induced 
pressures is not simple. Use of a typical acoustic filter would not 
be satisfactory since the filter connections would only furnish more 
nonuniformities with which to generate more noise. The use of a 
compliant section at the valve end of the supply transmission line 
might, however, prove to be useful. Also, low pressure systems are 
less susceptible to the noise due to their low density, and perhaps 
should be exclusively used where very high accelerations will be 
encountered. Securing the piping to stiff structure at short 
intervals would also reduce the transverse vibrations and make the 
system susceptible to vibrations in only one axis, longitudinal. 
Lastly, a thorough vibration test of the empty piping system should 
be made since high induced pressures can be expected at any frequency 
where strong structural resonances exist. 
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Figure 2.1 Basic Fluid System 



FLAPPER- NOZZLE 

RESERVOIR CONTROL VALVE 
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Figure 2.3 Pipe Element 
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Figure 2 4 Sample Problem 
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Figure 2*5 Pressure Frequency Response For Sample Problem 
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Figure 2.6 Shape of Damped Reaonant Feaks 
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Figure 2.7 Test Setup for Longitudinal Pipe Vibrations 




Figure 2.8 Test Setup for Transverse Pipe Vibrations 
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Figure 2„9 Experimental Arrangement 
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Figure 2»10 Pressure Frequency Response in a Pipe Vibrating Longitudinally 
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Figure 2.11 Velocity of Sound Frequency Dependence 
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Figure 2.12 Acceleration Frequency Response For a Pipe Vibrating 
Longitudinally 
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FIGURE 2.13 PRESSURE SPECTRUM IN A PIPE FORCED LONGITUDINALLY WITH 
WHITE NOISE EXCITATION 
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Figure 2^15 Explanation and Scale for Reading db 



db “ 10 log ^ where q is the reference quantity, 
^o 

however, are in units of quantity squared. Therefore, 



Spectrum, 



2 

db ~ 10 log ^ — 2 ~ 

q 

‘o 

Graphs that are presented in db can be easily read by using the 
scale below. The correct quantity for a given db is obtained by 
multiplying the reference quantity by the appropriate factor in the 
scale. Fcr units between decades, the decade factor is first multi- 
plied by the appropriate unit factor and then times the reference 
quantity, i,e, 23 db = 10x1.41 = 14,1 q^. 
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FIGURE 2.16 VOLTAGE SPECTRUM OF WHITE NOISE GENERATOR 
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FIGURE 2 23 PRESSURE SPECTRUM IN A PIPE FORCED BOTH LONGITUDINALLY AND 
TRANSVERSELY WITH WHITE NOISE EXCITATION 
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li. CONTROL VALVE ANALYSIS 

3.1 I NTRODUCTION 

In Chapter 2, fluid noise introduction through the environmental 
vibrations of a passive transmission line element was developed. In 
this Chapter an active element, a flapper -nozzle control valve, will 
be investigated while in the same vibrating environment. There are 
two transmission lines associated with the control valve, the supply 
line leading to the upstream orifice and the load line that transmits 
the signal to the load orifice. Noise can be introduced into the 
system in either or both of these. Induced relative motion of the valve 
flapper also creates a noise problem and will be investigated along 
with the two other sources. The small dimensions of a control valve 
prohibit the generation of structural induced internal pressures since 
the resonant frequencies are so high as to make the pressures 
insignificant . 

A quasi-linear approach will be used whenever nonlinearities are 
encountered by linearizing the relation around the operating point. 

This method gives satisfactory results for small excursions of the 
variables . 

In the following sections the valve chamber and flapper will 
be analysed first. Then two different models, distributed and lumped, 
will be constructed for the load line. Finally, the experimental 
results will be presented and discussed. 
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3.2 MATHEMATICAL MODELS 
3o2.1 VALVE CHAMBER 



In a flapper-nozzle control valve the variable orifice is the 
curtain area between the nozzle and flapper. This orifice is located 
downstream of a fixed orifice so that any change in its area alters 
the mass flow rate out of the valve and causes the pressure in the 
chamber to vary, see Fig. 3.1. For the valve analysed, both the fixed 
and variable orifices have choked flow through them and behave as pure 
resistances . 

Writing the continuity equation for the valve chamber of Fig. 3.1 

gives 



w 

n 



+ w. 



w 

s 




P V 
c c 



0 , 



(3.2-1) 



, .532N„A_ps j .532N_uD_x D„ u i j 

where w = ^ and w = n n tv-S. are the choked mass 

flow rates through the upstream supply orifice and the variable nozzle 

orifice respectively. In these relations the N’s are the orifice 

discharge coefficients, T is the absolute temperature, and A and 

o s 

ttD X are the areas of the respective orifices. w„ is the mass flow 
nr ^ I 

rate into the transmission line, and the valve chamber volume, V^, is 



taken to include that of the bellows. 



Let w. = 



.532N A p . 

s S SI 



SI 






R 



be the initial supply steady 



state mass flow rate at the design operating point, where R^ is the 



resistance of the orifice, and the subscript i is used to denote 
initial steady state condition. Then using the equation of state, 
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, and the differential operator, D, Eq . (3.2-1) can be 



R T 
o o 

nondimensionalized to give 

X p R w 
r^c , s Z 

^ ^ 

X .p 0 p . ^ . 

ri^si ^ SI ^si 

and for the steady state, 



5 p p . Dx 

® + tD — + Tb — = 0 (3.2-2) 
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SI ri 
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SI 



In these equations 
N 



N ttD X , 
n n ri 



N A 
s s 



V .p . 

Cl^Sl 

" = VTT- ’ 

1 o o 



A,x „p , 1 . 

b ri^si b 

’'b ~ w.R T a 
1 o o n 



where R is the gas constant and A. is the area of the bellows. The 
o b 

lengths Z. and Z are defined in Fig. 3.1. 
b n 

In order to get rid of the nonlinearity in the first term, 

Eq . (3,2-2) can be put into a quasi-linear form. Thus, 



p . Ax Ap AWp 
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SI 



tD -f N 



(3.2-4) 



In this equation it is seen that the valve chamber gives a first 
order lag to the pressure response. For optimum operation, the time 
constant t should be small and the signal to be transmitted introduced 
through the flapper displacement, Ax^. However, with the valve in a 
vibrating environment, noise can be introduced by the transmission lines 
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in the form of Ap^ and and also as an additional induced 

flapper displacement. 

Ap^ is the variation in pressure at the upstream orifice induced 

by the supply line vibrations. Since the orifice is choked, its 

impedence as seen by the line is a constant resistance and independent 

of the valve characteristics. Therefore, Ap can be independently 

s 

determined by the methods of Chapter 2 and applied as a forcing 
function in Eq. (3.2-4). 

On the other hand the flapper motion, Ax^ , is very dependent 
on the rest of the control valve variables. This is also true of 
For although it is a transmission line induced mass flow, it cannot 
be independently determined since the impedence as seen by the line 
is the valve chamber. This complicated impedence is wholly dependent 
on the other variables. Therefore, to completely describe the system, 
relations for the flapper and load line will have to be developed. 

3.2.2 VALVE FLAPPER 

Consider now the behavior of the flapper. The force of the air 
jet acting on it can be found by applying the noninertial momentum 



equation to the control volume in Fig. 3.2, and taking into account 



The noninertial momentum equation is 




where the r stands for relative to the accelerating control volume. 



the unsteady terms which previously have 




. This gives 
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''C.V, 

(3.2-5) 



where is the force of the flapper on the fluid (or the fluid on 

the flapper), m is the mass of air in the control volume, and x is 
a o 

the motion of the platform. 

The above equation was obtained by noticing that the pressure 
forces cancel each other ever3nA7here except over the nozzle area, A^. 
Also in the first term on the right hand side, it is assumed that the 
air flowing outward through the curtain area has a mean axial relative 
velocity of one half the flapper relative velocity. In the last term, 
the integrand vanishes everywhere except in the volume between the 
flapper and nozzle. In order to evaluate this term a simple approxi- 
mate, one-dimensional flow solution can be integrated over this volume* 
This has been done in Appendix B with the results that 
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Also from the equation of state, 

P v2 

2 Pn n , ^ 2 

rr = 



P V 
n n 



(3.2-7) 



o o 



where k is the ratio of specific heats, and is the Mach number. 

Since the chamber area is much larger than that of the nozzle, the 

relation for stagnation pressure can be used to give 
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^ . (1 + nh 

p 2 n 



2 



(3.2-8) 
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Thus substituting Eqs . (3,2-6), (3.2-7), and (3.2-8) into (3.2-5), 
one obtains 
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(3.2-9) 



The term in brackets is merely a function of k and and thus a 
* * 

function of A , where A is the curtain area for choked flow. 



n 



However, 



n 



^ 4x 
n _ I 

A *“ D 
n n 



Therefore this term can be replaced with an emperical expression that 
4x 



is a function of 



This was done by fitting a second order 



n 



relation to it, Fig. 3.3. The result, 
4x 



,243 ( -^ )^ + 1 
n 



4x 



gives excellent agreement up to — = .7. Virtually all flapper- 

n 

nozzle valve operation is well below this point. 

Now writing Eq. (3.2-9) in quasi-llnear form after substituting 
in the above relation, the differential operator, and w^, gives 
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the expression for the force the fluid exerts on the flapper, where 
4x 

( ^ 

X - ( D ^ 

n 

The flapper is positioned by a torque motor and bellows feedback, 
Fig. 3.1. Because the torque motor mass, damping, and stiffness are 
considerably greater than those of the flapper, a lumped model for the 
flapper can be used giving the system a single degree of freedom. 

From the Figure it can be seen that an inertial force will act on the 
flapper if the platform experiences accelerations in the direction of 
the control valve principle axis. This will also be true for the torque 
motor rotor if it is not S3rmmetrically balanced. 

The differential equation for the balanced condition is 
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(3.2-11) 



where F^ is the signal torque applied to the torque motor and and 
k^ are the stiffnesses of the bellows and spring respectively, m^, 
j^, b^, b^, k^, and k^ are the mass, inertia, damping, and stiffness 
of the flapper and torque motor, and the lengths ly and 

are defined in Fig. 3.1. When this equation is written in the quasi- 
linear form, put into Eq. (3.2-10) and the whole expression nondimen- 



sionalized by kx . , one obtains, 
^ ri 



AF - AI = 



11 



-2 + (2.; + =1 ns 

03 ^S1 

n 



ci V D 



) — +(1 + K - .485aX 



ci. 



SI 



Ax 



ri 



■7 nN — - a( .243X + 7^+1) 
3 0) A 









(3.2-12) 



M I 




72 



In this Equation, 
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Eq , 0*2-12) shows that the flapper motion is dependent on the 
chamber pressure. It also shows that the pressure is a function of 
the inertia force and susceptible to noise introduction by it. 

Still one more equation is needed in addition to (3.2-3), (3.2-4), 
and (3.2-12) to completely describe the system. This is the relation 
for the load transmission lineo Models for this line will be constructed 
in the next two sections. 



3.2.3 TRANSMISSION LINE , DISTRIBUTED PARAMETERS 

The transmission line leads from the valve chamber to the choked 

load orifice. Depending on the vibration field, it can be excited 

longitudinally, transversely, or in a combination. Regardless of the 

excitation or the piping supports, the methods of Chapter 2 can be 

used to obtain a relation for Aw^ in terms of Ap^. Here the pipe will 

be assumed rigid and forced, along with the valve, in a longitudinal 

2 

direction with an acceleration of D Ay . 

•^o 

Using the relative transmission matrix from Section 2.4 for the 



fluid, one obtains 
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Ap = cosh r Ap(L) -f Z sinh T Aw (L) + Z p ADsinh T Ay 
^ c r CO “^o 



and 



(3.2-13) 



AWj^ = — sinh F Ap(L) + cosh F Aw^(L) + p^AD(coshF- l)Ay^, 
c 

where Aw^ is the relative mass flow rate out of the chamber, Eq , (3.2-1). 

Since the load orifice is choked, w (L) = ^ 4' ^ ' > where R is the 

r K 

orifice resistance. Therefore, 
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(3.2-14) 



and in the steady state. 
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Putting this last expression into Eq. (3.2-3) gives. 
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(3.2-15) 



Eqs. (3,2-14) and (3.2-15) along with (3.2-4) and (3.2-12) 
completely define the system. Eq. (3.2-14) is for one specific type 
of load line and excitation while the other three Equations are general. 
The system behavior will be analysed using these Equations in Section 3.3, 
Before going on to a lumped model development for the transmission 
line, it is of interest to consider the impedence the line sees looking 
into the valve. Once this is known the load transmission line is 
completely determined and can be independently analysed. The impedence 
only has meaning when the valve is passive or AF, AI, and Ap^ are zero. 



For this case, since Z = - — 

v Aw 

I 

and (3.2-15) gives 
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, combining Eqs. (3.2-4), (3.2-12), 
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From Eq. (3.2-16), the impedence and real and imaginary parts are 
plotted in Fig, 3.4 for the particular control valve used in the 
experimental work. Also included is the impedence for when there is 
no bellows feedback, or = k = == 0. It is seen that feedback 

emphasizes the flapper resonance, but that it decreases the impedence 
at lower frequencies. 

It is to be noted that the real part of the impedence is positive 
up to the resonant frequency where it becomes negative, and the 
imaginary part is negative for all frequencies. Using this figure 
along with the known orifice resistance, Eqs. (3.2-13) can be used to 
solve the load line pressure and mass flow rate at any frequency. 



3.2.4 TRANSMISSION LINE, LUMPED PARAMETERS , AND ANALOG COMPUTER 
MODELLING TECHNIQUES 



In addition to the distributed model for the line, it is also 
desirous to construct a lumped model which allows the system to be 
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analysed on the analog computer. The advantage of using the computer 
is that each of the parameters can be easily changed and its effect 
on the rest of the system readily seen. 

The line is modelled as lumped inertias, resistances, and 
compliances. The inertias and resistances are combined into series 
impedence elements and the compliances constitute shunt admittance 
elements, see Fig. 3.5. 

Causality must be taken into account when constructing the model 
so that the problem can be solved on the analog computer without using 
differentiation which causes difficulties. The specified variable, 
pressure or mass flow rate, at each end of the line determines what 
the nearest lumped element must be. Pressure requires a series 
impedence while mass flow rate requires a shunt admittance. Also, the 
number of line resonances a person desires to model determines the 
number of inertia elements that must be used. Thus since the trans- 
mission line in this investigation starts with a known pressure, Ap^, 
and terminates with an induced flow source due to the motion of the 
line, and since it is desired to model the first two resonances of the 
line, the model in Fig. 3.5 is appropriate. 

The transmission matrix for the first half of the line is 
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and that for the second half is 
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respectively, when the substitution Aw (L) = is made. 

r K 

Now extending the results in Section 2.4, the series impedence 
becomes 

goA £ 

= Z — + R 
c c £ 



= Z T D + R, 
c e I 



(3.2-17) 



where T = ^ is the transmission time for a wave to travel the 
e c 

length of the pipe, and R. is the total frictional resistance. Also, 



the shunt admittance becomes 
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(3,2-18) 
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Putting Eqs, (3.2-17) and (3.2-18) into the transmission matrices 
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For the steady state, Eqs. (3.2-19) give 
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which when combined witli Fq . (3.2-3) becomes, 
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(3.2-20) 
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Although Eqs. (3.2-19) could be cascaded to give one overall 
transmission matrix, this should not be done for the same reason as 
mentioned previously, namely the resulting equations would not be 
able to be solved on the analog computer without the use of differenta- 
tion. Therefore, regardless of the number of lumped sections the line 
is divided into, each section should be handled separately. 



The resistance is not included in the forcing term since only 
the flow relative to the pipe causes a loss of pressure. 
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In order to determine for use in the above equations, 
necessary to use a quasi-linear approach since the frictional 
tance is nonlinear. The pressure loss equation is 
4fLp 

o 



2g D, 

where f is the friction factor, p is the mean fluid density. 



is the pipe diameter. Now putting in mass flow rate gives 
2 



P = 



4fLp w 
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where A is the flow area. 



Putting this into the quasi-linear form it becomes 



Ap = 



SfLp^w^Aw 
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( 3 . 2 - 21 ) 



Using Eqs. (3.2-4), (3.2-12), (3.2-19), and (3.2-20) it is now 
possible to set up the problem on the analog computer. An operational 
block diagram for this setup is contained in Fig. 3.6. Since every- 
thing is nondimensionalized , including time, scaling is not too diffi- 
cult and is indicated on the Figure. The various inputs are also 
marked as well as the system variables that are of interest. It is 
noted that the load line induced noise input, enters in at two 

different points due to the load line being modelled in two parts. 

This analog model now allows one to vary the input and system 
parameters, and to investigate their effect on the system variables. 



Results thus obtained are discussed along with the other analytic 
and experimental results in the following section. 



3.3 EXPERIMENTAL WORK AND DISCUSSION OT RESULTS 

A small flapper-nozzle control valve was built in order to 
verify experimentally the analytical results. It was equipped with 
a removable bellows feedback system, Fig. 3.7. The flapper displace- 
ment could be measured by means of a capacitance pickup, and the 
chamber pressure and supply line pressure just upstream of the ori- 
fice were obtained with a piezoelectric transducer. The flapper 
displacement was controlled by means of a torque motor supplied 
from a servo amplifier. A view of the general arrangement and instru- 
mentation are shown in Fig. 3.8, and further discussion is contained 
in Appendix C. 

Since a vibration generator large enough to shake the whole 
system was not available, it was necessary to simulate the noise 
inputs. The pressures in the supply line were generated by supplying 
the valve through the load line of another flapper-nozzle valve 
operating at a higher pressure. In this manner pressures at any 
frequency could be imposed on top of the mean supply pressure of 
100 psia, simulating resonant noise peaks. 

The behavior of the control valve transfer function of, 
when subjected to supply line noise, is shown in the frequency 

response curves of Fig, 3.9. The Figure contains an experimentally 
obtained curve for the valve operating with bellows feedback as well 
as theoretical curves with and without feedback. 
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The theoretical results were obtained by combining Eqs. (3.2-4), 
(3.2-12), (3.2-14), and (3.2-15). After putting in the appropriate 
values for the valve used, this distributed model equation becomes 
aAp + .59bAG. + .012acAG 

Ap = TT^-T (3.3-1) 

c ae + 6.75b + 1.82ad 



where Ap^ is in psi, Ap^ is the supply line pressure noise in psi, 

AG^ is the acceleration in g’s experiences by the valve in the direction 
of its principle axis, and AG^ is the acceleration of the load line 
in g’s. Any one of the three noise transfer functions can be ob- 
tained by setting the other two noise inputs to zero. In this 
equation, 
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When there is no bellow feedback, the denominator of Eq. (3.3-1) 



becomes ae + 1.30b + 1.82ad, 
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and a = .35 — i + 1.00 - — , 
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b = .22, 
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In the Figure, the transfer function is seen to behave as a 
first order lag cascaded with a second order system. This is due to 
the valve chamber and flapper respectively* Feedback increases the 
effect of the flapper resonance, but it decreases the pressure at 
lower frequencies. The experimental curve is not continued past 
700 cps because the response of the first (noise producing) flapper - 
nozzle valve drops off to zero at this frequency, and no noise could 
be introduced to the test valve. 

The system appears to have slightly more damping than was 
estimated, and therefore the experimental peak is not as high as 
predicted, and it is at a slightly lower frequency. From the Figure, 
the damping is calculated to be approximately .08 and .10 for the 
theoretical and experimental results respectively. It is to be noted 
that the flow noise in the valve created a region of uncertainty in 
the experimental results having a width of approximately .02 on the 
Figure . 

It appears that if the flapper-torque motor arrangement had less 
damping and if the bellows feedback were increased, then t^e noise 
transfer function would become much larger near the flapper resonance. 
In order to study the effect of varying these two parameters, analog 
computer solutions were obtained using the lumped model formulated in 
Section 3.2.4. 

It was found that the noise signal could be increased as expected. 
However of much more interest, it was also found that the total noise 
power could be amplified in the valve at frequencies near the flapper 
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resonance. That is to say that the noise power associated with the 
bellows, force times velocity, was greater than the noise power 
entering the valve, pressure times flow. This amplification can take 
place because some of the steady power flowing into the valve is 
converted into noise power. This is a danger of noise in an active 
element and differs from what can happen in a passive element. In a 
passive element high noise signals can be obtained through resonances 
as in the transmissian lines studied, but the total noise power cannot 
increase. It is of interest that the increased power was measured in 
the bellows which is a direct analogy to a load piston operating with 
the control valve. Thus it is seen that the noise can be amplified 
into actually displacing a load. Also, increased motion of the 
flapper due to the feedback of the supply line noise can cause the 
flapper to strike the nozzle. This flapper motion saturation makes 
the valve inoperative to command signals. 

In terms of numbers, the threshold of amplification was observed 
when damping was decreased by a factor of .7 and feedback was increased 
by a factor of 1.7. Another combination was found when damping was 
decreased by a half and feedback was increased by a factor of 1.2. 

These factors decreased the flapper system damping ratios to .25 and 
.17 respectively for the two cases. These are not excessively low 
values and might well be encountered in many systems. These results 
could not be checked experimentally using the test valve because 
the feedback and torque motor damping had fixed values. When the torque 
motor was removed to decrease the amount of damping, the system became 
unstable. Since the feedback could not be increased, it also was not 
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possible to make tbe flapper motion saturate. The maximum flapper 
motion that could be excited was about .0015 inches while the steady 
state nozzle clearance was .005 inches. 

In order to simulate inertial excitation to the flapper, a white 
noise signal was added into tbe servo amplifier input. This showed 
the first resonant peak as expected, but it also showed smaller 
peaks at the equalizing spring resonances and at the second resonance 
of the flapper. Fig. 3.10 is the frequency response transfer function 
of the chamber pressure per g of induced flapper excitation. The 
characteristics are very much like those discussed in conjunction 
with Fig. 3.9. There is, however, a small peak at 300 cps in the 
experimental result due to the equalizing spring, and the pressures at 
the lower frequencies are larger than predicted. Increasing feedback 
and decreasing damping would have the same effects with regards to 
amplifying this induced noise as was discussed with supply line noise. 

It is also possible to get flapper motion saturation. 

Since there was no driver available that would work against the 
pressure in the load line, it was not possible to simulate load line 
induced vibrations. The line was hooked up to the vibration generator, 
but this did not induce signals large enough to be recorded. Fig. 3.11 
shows that the theory predicts that these pressures will be very small. 
Even for 100 g*s of excitation the induced pressure would be less than 
one psi even at the flapper resonant frequency. The smallness of the 
induced pressure is due to the capacitive effect of the valve chamber 
volume, anH to the low pressure of the load line, approximately 35 psia, 
causing the air density to be quite low. A valve operating at higher 
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pressures would be much more susceptible to this type of noise. 

The above results have shown that the control valve is susceptible 
to noise from the supply line and inertial forcing of the flapper. 

The noise amplitude is greatest at low frequencies and near the flapper 
resonant frequency. When the flapper system damping is low enough, 
bellows force feedback can cause amplification of the noise power 
to the valve and saturation of the flapper motion. 

External vibrations of the load line do cause resonant fre- 
quency peaks, but because of the reduced pressure of this line, it is 
only of concern in valves with very high supply pressures. 



3.4 PROCEDURE FOR TESTING CONTROL VALVE SUSCEPTIBILITY TO NOISE 



Of considerable interest to the person who must specify a parti- 
cular control valve for a given job, is whether or not the valve will 
be highly susceptible to vibration induced noise. Therefore, it should 
become standard procedure to test any such contr^^l valve before it is 
assigned for use. Since there are three ways noise can be induced 
into the system, there are primarily three tests. 

The first test deals with noise introduced through the supply 
line. For this test, the supply line noise can be generated in the 
same manner used in this theses . Supply the valve through the load 
line of a second control valve operating at a higher pressure. In this 
way sinusoidal noise signals can be imposed on the mean supply pressure 



for any frequency. The frequency response for \ t ^ I should then 

Ax ^s 

be obtained as well as the response for I > making sure to test 
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at the maximum feedback that will be used in practice. If either of 

these two responses is excessively large, the control valve should not 

be used where environmental vibrations could be a problem. 

The second test is for vibration induced flapper motion. If a 

large vibration shaker is available the valve can be vibrated directly 
Ax Ap 

and 1"^^ (co)|and|— ^ (lo) | recorded. Otherwise the frequency response 

for both of these can be obtained by applying a signal to the torque 

AF^ 

motor as was done in this thesis. In this case G = where F 

£m^ t 

is the torque applied to the motor, m^ is the equivalent forced mass 
of the flapper, and i is the lever arm of the rotor. If either of the 
responses is too large, an effort can be made to balance the rotor and 
flapper assembly. If one is fortunate enough to have the large shaker, 
any reduction in response due to balancing will show up directly. 
However, if one is simulating the vibration, balancing will reduce the 
value of m^ used in calculating AG. If m^ is reduced enough, AG will 
be excessive in order to force a response, and the control valve will 
be safe for use. If balancing is not successful the valve should not 
be used. 

The third test is for load line induced pressures. If the control 

valve chamber does not operate above several hundred psla, this type 

of noise should not be a problem. A test can be made with a small 

vibration generator used to vibrate the end of a flexible line in a 

Ap 

I * C I 

(co) I , recorded. The 

amplitude of this response will dictate whether the valve is safe or 
not . 

One additional test can be conducted with the fluid transmission 
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lines, A thorough vibration test of the empty piping system should 
be made since high induced pressures can be expected at any frequency 
where strong structural resonances exist. If any of these resonances 
are near the flapper resonant frequency of the proposed control valve, 
a different valve should be used. 






• # 



I 




4 











87 



Figure 3.1 Control Valve Arrangement 
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Figure 3.2 Flapper-Nozzle Control Volume 
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Figure 3.3 Nozzle Relation Approximation 




Figure 3.4 Valve Chamber Impedence as Seen by the Load Line 
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Figure 3.5 Lumped Parameter Model of Load Transmission Line 
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Figure 3.6 Analog Computer Operational Block Diagram 
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Figure 3.7 Flapper-Nozzle Control Valve 




Figure 3.8 Experimental Setup 




Figure 3.9 Chamber I’ressure Frequency Response to Supply Line Noise 
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Figure 3.10 Chamber Pressure Frequency Response to Induced Flapper Motion 
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Pigur* 3.11 Chamber Fressurc Frequency Response to Load Line Longitudinal 
Vibrations 
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4. CONCLUSIONS 

4.3 VIBRATION INDUCED TRANSMISSION LINE NOISE 

From the results of Chapter 2 it is concluded that piping 
subject to white noise environmental vibrations can excite resonant 
pressure peaks in the contained fluid over a wide band of frequencies. 
The pressures are caused by induced fluid flow at pipe boundary 
nonuniformities (bends and area changes) due to the relative motion 
of the piping. The peaks occur at the various fluid and piping 
resonant frequencies. 

Either longitudinal or transverse piping support vibrations, 
or a combination, will cause induced pressures. In all cases it is 
the longitudinal motion at the end nonuniformities that forces the 
fluid. This is true even for transverse vibrations since the trans- 
verse pipe motion excites longitudinal vibrations of the ends at' the 
pipe bending resonant frequencies and at double the resonant 
frequencies. Since pressures are induced for either type of vibra- 
tion, a fluid piping system is susceptible to noise generation 
regardless of its orientation to the vibration field. 

The level of induced pressure is proportional to the fluid 
density, the effective speed of sound in the fluid, and the amount 
of excitation. Resonant peak levels of .6 psi per g of longi- 
tudinal acceleration can be obtained at the lower frequencies for 
air at 200 psia mean pressure in one quarter inch copper tubing. 
Preliminary results show that correspondingly higher pressures can 
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be obtained with liquids due to their higher density* The 
pressures generated by the first several pipe resonant frequencies 
when subjected to transverse vibrations are comparable to the 
longitudinally excited case. However, levels due to transverse 
double resonant frequency excitation are less by about a factor 
of two . 

Extremely high pressures can be excited when a fluid resonance 
coincides with that of the piping system. For a more complicated 
system it is impossible to distinguish between these two types of 
resonances, but there are still corresponding absolute system 
resonances th^t are higher than the others. 

4.2 CONTROL VALVE 

:fe 

For the one type of active element studied while in a 
vibrating environment, noise introduction into the valve takes 
place through the supply transmission line, induced motion of the 
flapper assembly, and for high pressure systems, through vibrations 
of the transmission line to the load. Since the valve transfer 
characteristic is that of a first order lag cascaded with a second 
order system, due to the capacitive chamber and flapper resonance 
respectively, the noise transfer functions are largest at the low 
frequencies up to the break frequency and again near the flapper 
resonant frequency. 

* 

A Pneumatic nozzle-flapper control valve. 



* 
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Force feedback in the way of a bellows arrangement reduces the 
pressure levels at the lower frequencies. However, it also intensi- 
fies those near the flapper resonance. By Increasing the feedback 
enough, it is possible to get amplification of the noise power for 
frequencies near the flapper resonance. That is to say that in any 
active element it is possible to convert some of the steady flow 
power entering the valve into noise power. This phenomenon does 
not occur in a passive element where signal amplification can take 
place, but not power amplification. 
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Ll, recommendations 

5.1 FUTURE WORK 

Since the results of this study have shown that significant 
vibration induced pressures can be generated in a fluid control 
system, further work in this field is warranted. The most fruitful 
areas for future work are the following. First, it is recommended 
that extensive experimental work with transversely vibrating pipes 
be conducted in order to better understand the structural coupling 
between longitudinal and transverse motion. Data should be taken 
for various pipe end constraints and different support geometries. 
Secondly, different types of active elements should be analysed 
while operating in a vibrating environment. This should include pure 
fluid elements such as bistable and Jet amplifiers, as well as systems 
with several stages of feedback and amplification. The results should 
be compared along with those for the flapper-nozzle valve to determine 
which elements are least susceptible to vibration induced noise . 

Further experimental work should also be conducted with the flapper- 
nozzle valve to verify that noise power amplification takes place. 
Lastly, it is recommended that parallel work be conducted for a 
hydraulic system. Preliminary work has shown that a liquid system 
will be even more susceptible to vibration induced noises. 

5.2 APPLICATION 

Information learned in this thesis can be applied in order to 
reduce a system's susceptibility to noise. Thus, it is recommended that 
all fluid equipment that will be exposed to strong environmental 
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vibrations be thoroughly tested before use. The system should be 
shock mounted if possible. Otherwise it should be mounted on stiff 
supports so that all resonant frequencies will be high. The 
transmission lines should also be secured at short intervals in order 
to decrease transverse vibrations. The installed system, empty 
of fluid, should then be vibration tested. High induced pressures 
can be expected at any frequency where strong structural resonances 
exist. Lastly, the fluid amplifiers should be carefully tested in 
order to pick those with good transmission characteristics that will 
not cause noise amplification. 
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APPENDIX A 

TRANSMISSION MATRICES 
A.l TRANSMISSION MATRIX DEVELOPMENT 



* 

A model for a uniform, one-dimensional, distributed system is, 




where Z is the series impedance per unit length, and Y is the shunt 
admittance. Either or both can contain the operation D = . 

Using pressure, p, and mass flow rate, w, as variables (force and 
displacement for pipes) the relations for a micro-element become. 



or 



dx 



- Zw 



^ - Y,^ 

dx " 



dx 



p(x) 




0 


Z 




p(x) 


w(x) 




Y 


0 




_w(x) 



(A. 1-1) 



References [isj and give extensive developments containing the 

material covered in section A.l, 

When working with mechanical and fluid systems, this form of the 
equations is obtained when the pressure or force variable is taken as 
compressive and the mass flow rate or displacement variable is taken 
in the direction of x, i,e. — ^x 
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Ordinary differential equations are written instead of partial 
because time is contained in the operator D in the system parameters* 
Writing the above matrix equation in compressed form and using 
matrix algebra, one finds that 



-ASw 

dx 

-Ax 

S(x2>=e 

-Ax, 

S(x^) = e 



and 



/T '^''2 C 
C = e S (x J . 



Therefore 






(x^) 



where ([^ is the integration constant matrix. 

* 

The transmission matrix can be expanded into, 



A(^2" 



(A. 1-2) 



cosh y(x 2 ~ Xj^) 



sinh y(x 2 ~ x^^) 



~ sinh y(^ 2 ~ y(x 2 - x^^) 



where Z^= is defined as the characteristic impedence and 



Y = as the propagation operator. 

The final expression becomes. 







cosh y(x2~ x^) sinh 




P(X2> 


w(Xj^) 


s: 


^ sinh y(x^- x^^) cosh rCx^- x^) 
_ c 




w(x2> 



(A. 1-3) 



This can be done by using either the Cayley Hamilton or Sylvester 
theorems [2^ . 
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This transmission matrix is applicable to any uniform, one-dimensional, 
linear, distributed system of two variables. When working with piping, 
p and w are replaced by f, force, and £, displacement, respectively. 



A. 2 TRANSMISSION MATRIX USING RELATIVE FLOW OR DISPLACEMENT VARIABLES 



There are times when working with a system mounted on a vibrating 
platform when it is more convenient to use flow or displacement variables 
relative to the platform motion. Thus for distributed fluid and 
mechanical systems the relations have the form 



= - Zw - Zw 
dx r o 



dw 
I 

dx 



- Yp 



(A. 2-1) 



since w(x) = w (x) + w . The r and o stand for relative to the 
r o 

platform and platform respectively. The platform is rigid so that 
its motion is independent of position. 

Using matrix algebra one obtains from 



that 



dx 



p(x) 




0 Z 




p(x) 




- Zw 

o 




= - 








+ 




W^(x) 




Y 0 




w (x) 
r 




0 



dS(x) 

dx 



= - (x) +([] 



S(xJ = 



-/^X„ -1 

- 'Co-A\C 



S(x^) = e 



e 



1 




I 
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Therefore 




A 



C^2~ 



Xf) 



S (^2) 




(A. 2-2) 



Now since 

a" 



0 1/Y 




- Zw 

rs 




and ^ = 


\J 


1/Z 0 




0 



the final expression becomes, 



p(x^) 




cosh y(x 2 “Xj^) 


Z^sinh y(x 2 ~Xj^) 




P(x2> 






1/Z^sinh y(x 2 ~Xj^) 


cosh y(x 2 ~Xj^) 




(x2> 



•f 



Z sinh y(Xo“X-)w 
c 2 1 o 



(A. 2-3) 



j^cosh y(x 2 ^x^)- 1 J 

It is to be noted that Eq, (A. 2-3) can be obtained more directly by 
substituting the expression w(x) = w^(x) + into Eq, (A.1-3) • This 
method is not recommended, however, for more complicated systems 
involving matrices of higher order. In those cases, application of 
Eq. (A, 2-2) will be more profitable. 



A. 3 TRANSMISSION MATRIX WITH TRANSPORT 



:k 

By introducing wave-scattering variables one can find a good 

approximate form of the transmission matrix when there is transport 

*For a comprehensive presentation of wave-scattering variables see 
Paynter |^1^ , pp. 268-271, and Brown . 
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present. This would be useful for cases of mean flow in fluid lines. 
Wave-scattering variables can be obtained from the previous expression 

d S 

Z — ” obtaining the canonical form of the differential 

^ -1 
equation by diagonalizing /\ . Thus, premultiplying by gives 



--n'^AAA''s. 

Defining S as S = ^ S 



-1 



leaves /^~ ^ as the desired diagonal matrix. 

Therefore, for a two variable system, 



and 



A 



n 



-1 1 



rr 



c c 

I'VT w 

C ’ C 



-■ 






u(x^) 



v(xi^ 






,-y(x2-Xi) 



u(x2> 



v(x2> 



(A. 3-1) 



Since 






and 



u(x) 

v(x) 



-1 



= S (x) = S (x) . 
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For a wavelike transmission line, i.e. fluid line, Eq. (A. 3-1) 
becomes , 



u(x^) 




v(x^) 





0 



-TD 



u(x^) 



vCx^) 



(A. 3-2) 



where T is the wave propagation time between points 1 and 2 and is 

X2-X1 

equal to — - — where c is some characteristic velocity of the wave 
propagation. When one looks closer it is seen that the scattering 
matrix contains only a pure time delay and a pure advance operator. 
Thus writing Eq. (A. 3-2) in the time domain it becomes, 

U2(t) = Uj^(t-T) 

v^(t) = v^Ct-T)^ 

This means that u is a rightward traveling wave and v is a leftward 
traveling wave 






A 



v^ 



as u^ is u^ delayed by T, and v^^ is ^2 delayed by T. 

Now if one introduces a mean transport velocity V, 



V u. 



the propagation times will be different for the two directions. For 

X2-X1 

the rightward traveling wave, ^ ^ , and for the leftward traveling 

X -X 



Now Eq . (A. 3-2) becomes 
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and transforming back to the S variables it becomes, 
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VT' VF" 
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T D 
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w(Xj^) 
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^ c c 
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-T^D 
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_ C 


-YT 
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Multiplying gives 







T D -T.D T D -T.D 




— - 


P(Xj^) 




r. Z — ^ 1 C 

e+e Z(e-e ) 

c 




r(x2> 




_ 1, 
0 


T D -T„D T D - T„D 






w(Xj^) 


z 


L {e. - e ) e +e 




w(x^) 



Now one can make the approximations that 



X2-X1 x^-x^ 



c + V 

c(l+-) 

c 



(1 _ ^ ) 
c c 



and 



T. = 



X2 



I c - V V . 

c(l - - ) 



1 (1 - ^ ) 
c c 



since — << 1 for all cases where the above analysis holds, 
c 



Substituting these relations back into the expression gives, 



P(xj^) 


-V(x2-Xj^)D 


cosh y(x 2 ~Xj^) 


Z sinh y(Xo-x-) 
c Z 1 




P(x2> 


w(Xj^) 


c2 


1/Z^sinh y(x 2 ~Xj^) 


cosh 




w(x2) 



(A. 3-3) 
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Thus it can be seen that transport merely introduces a multi- 
plicative factor and that Eq . (A. 3-3) approaches Eq • (A. 1-3) as 
V — ^0. 

When the equations are written in terms of a relative flow or 
displacement variables one obtains. 





-Vfx^-x^D r 


cosh y(x 2 “Xj^) 


Z^sinh 




W (x- ) 

V y 


1 


1/Z^sinh y(x 2 ~Xj^) 


cosh y(x 2 ~Xj^) 





+ 



Z sinh y(x^-x-)w 
c Z L o 

V(x2-x^)D 

P- 

jzosh yCx^^Xj^)- e 




(A. 3-4) 



In this case one of the internal terms is altered slightly. Yet this 
also approaches its corresponding no transport equation, Eq. (A. 2-3), 



as V— ^0. 
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APPENDIX B 

UNSTEADY MOMENTUM TERM EVALUATION 
B.l APPROXIMATE SOLUTION 

A simple one-dimensional flow solution can be integrated over 
the volume between the nozzle and flapper in order to obtain an 
approximate value for the integral. 



This is done by assuming that the one-dimensional flow lines are 
always perpendicular to the slant areas of the right circular and 
truncated cones into which the volume can be divided, Fig. B.l. 

The integration of this volume is done in two parts. The first 
includes the circular cone volume elements starting at the nozzle 
opening, and the second consists of the truncated cone elements 
which end at the curtain area. Dividing the problem this way gives 
an exact integration of the volume for the first part and a close 
approximation for the second. 

Beginning with the cones, let A^ = ttsR be their slant areas 
where R is the radius and s is the slant distance. Then an element 
of volume is 




A sd0 
s 



where -j is the factor for a cone, and 0 is defined in the Figure. 



It is also seen that 



I 
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R 

cos 6 
w 



Thus, since p = 



V A 
r s 



, these expressions can be put into the first 



part of the integral to give 




V COS0 w A R d0 
r r s 

g V A 3 COS0 
^o r s 
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— . Integrating then gives 
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For the truncated cone elements, A^ « 7Ts(R+r) is the slant area 
and r is the radius defined in the Figure. An element of volume is 
then once again 



dV = 



A sd0 
s 



where s = 



sin0 

Thus the remaining part of the integral becomes, 

TT 
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2 V COS0 w A X 
r r s r 

g V A sin0 
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9_ ('2 Vr 

3t 



de . 



Integrating, one obtains. 
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Since — is usually about .1 to .3, this expression can be approxi- 
K 



1 . W X X 

mated by t ), being Jin — - 
d 9t g^ 



- 1. 



Therefore, adding the two parts of the integral together gives 
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xr 3 



3t 
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since the relative flow through any slant area is equal to the flow 



into the nozzle, or w = w • 
^ r n 
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Figure B.l Elemental Volumes 
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APPENDIX C 

COMMENTS ON EXPERIMENTAL PROCEDURES AND INSTRUMENTATION 
C,1 TRANSMISSION LINE NOISE EXPERIMENTATION 

In Chapter 2, the physical set up for experimentally determining 
the induced pressures in vibrating pipes was described and figures 
provided. This Appendix contains some additional comments. 

While taking data it was necessary to secure the pressure trans- 
ducer to a nonvibrating support due to the sensitivity of the trans- 
ducer to vibrations. This was accomplished by having the copper tube 
terminate in a flexible section of poly-flow tubing which could then 
be secured. The other end of the piping was attached to the constant 
pressure reservoir in the same manner. 

The vibration generator was made from a 35 watt speaker unit 
with the cone removed. Extra care was necessary with alignment in 
order to keep the coil from rubbing while it was in operation. This 
speaker unit made a very satisfactory vibration shaker that responded 
well to input signals from both a sinusoidal oscillator and a random 
noise generator. These input signals were preamplified by a regular 
high fidelity amplifier before being fed to the shaker. It was 
necessary to monitor the current going to the shaker coil in order 
to keep from burning it out. 

The mean test pressure level was maintained by use of a regulator 
reducer that was fed from a high pressure air line. Any fluctuations 
were damped out by the constant pressure reservoir. A globe valve was 
installed upstream of the reservoir in order to bleed down the system. 
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The pressures and accelerations were picked up with Kistler, 

Model 601-A and Endevco, Model 2215 piezoelectric transducers 
respectively. The signals then went through buffer amplifiers to the 
oscilloscope and analyser and recorder. The frequency analyser and 
level recorder were Briiel and Kjaer, Type 2107 and Type 2305 
respectively. This equipment automatically swept through the frequency 
range with a 6% filter bandwidth and recorded the mean square levels. 
Averaging time was controlled by the paper feed rate in the recorder. 



C.2 CONTROL VALVE NOISE EXPERIMENTATION 

The flapper-nozzle valve built for experimental purposes had a 
chamber volume of .1 cubic inches. The nozzle diameter was .116 
inches while those of the supply and load orifices were »086 and .116 
Inches respectively. The valve was designed to operate at 100 psia 
mean supply pressure and with a steady state flapper displacement of 
.005 inches. 

The flapper displacement was controlled by a Midwestern, Type 
11 F torque motor that operated from a 300 volt d.c. American Measure- 
ment and Control servo amplifier. Model 652. The flapper displacement 
was measured by using a Decker Unit, Type 902-1. This capacitive 
pickup had the advantage of not imposing any unwanted force on the 
flapper . 

Pressure measurements were obtained with the same piezoelectric 
transducers previously mentioned. Readings could be taken in the valve 
chamber and just upstream of the supply orifice. The supply air was 
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provided through another flapper-nozzle valve whose load was the 
test valve* In order to maintain this mean supply pressure at 
100 psia, it was necessary to operate the first valve supply pressure 
at 467 psia due to the large mass flow rate required by the test 
valve* 

Data analysis was accomplished with the same equipment mentioned 
in the previous section* 
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